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coupling function based on renormalization- group considerations. 
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I. INTRODUCTION 



It is observationally confirmed not only that inflation occurred in the early universe, 
but that the current expansion of the universe is accelerating [1, 2]. Although there exist 
various scenarios to account for the late-time acceleration in the expansion of the universe, 
the mechanism is not well established yet (for recent reviews, see [3-7]). 

The scenarios to explain the late-time acceleration of the universe fall into two broad 
categories [6]. One is general relativistic approaches, i.e., dark energy. The other is modified 
gravity approaches, i.e., dark gravity. As the most promising one of the latter approaches, 
the modifications to the Einstein-Hilbert action, e.g., the addition of an arbitrary function 
of the scalar curvature to it, have been studied (for a review, see [7]). Such a modified 
theory is considered as an alternative gravitational theory, so that it must pass cosmological 
bounds and solar system tests. 

Recently, Hu and Sawicki have proposed a very realistic modified gravitational theory 
that evade solar-system tests [8] (for related studies, see [9]). In this theory, an effective 
epoch described by the cold dark matter model with cosmological constant (ACDM), which 
explains high-precision observational data, is realized as in general relativity with cosmolog- 
ical constant (for a review of observational data confronted with modified gravity, see [10]). 
Although this theory is successful in explaining the late-time acceleration of the universe, 
the possibility of the realization of inflation has not been discussed in Ref. [8]. In Refs. [11- 
13], therefore, modified gravities in which both inflation and the late-time acceleration of 
the universe can be realized, following the previous inflation-acceleration unification pro- 
posal [14], have been presented and investigated. The classification of viable F(R) gravities 
maybe suggested too [12]. Here, F(R) is an arbitrary function of the scalar curvature R. 

As another gravitational source of inflation and the late-time acceleration of the universe, 
a coupling between the scalar curvature and matter Lagrangian has been studied [15, 16]. 
Such a coupling may be applied for the realization of the dynamical cancellation of cosmo- 
logical constant [17]. The criteria for the viability of such theories have been considered in 
Refs. [18-20]. Recently, simple coupling between the scalar curvature function 

and the kinetic term of a massless scalar field in viable modified gravity has been consid- 
ered [21]. 

On the other hand, it is known that the coupling between the scalar curvature and 
the Lagrangian of the electromagnetic field arises in curved spacetime due to one-loop 
vacuum-polarization effects in Quantum Electrodynamics (QED) [22]. Such a non-minimal 
gravitational coupling of the electromagnetic field breaks the conformal invariance of the 
electromagnetic field, so that electromagnetic quantum fluctuations can be generated at 
the inflationary stage even in the Friedmann-Robertson- Walker (FRW) spacetime, which 
is conformally flat [23-25]. They can appear as large-scale magnetic fields at the present 
time because their scale is made longer due to inflation. These large-scale magnetic fields 
can be the origin of the large-scale magnetic fields with the field strength 10~ 7 -10~ 6 G on 
lOkpc-lMpc scale observed in clusters of galaxies [26] (for reviews of cosmic magnetic fields, 
see [27]). 

In the present paper, we consider inflation and the late-time acceleration of the universe in 
non-minimal electromagnetism, in which the electromagnetic field couples to the function of 
scalar curvature. We show that power-law inflation can be realized due to the non-minimal 
gravitational coupling of the electromagnetic field, and that large-scale magnetic fields can 
be generated due to the breaking of the conformal invariance of the electromagnetic field 
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through its non-minimal gravitational coupling 1 . The mechanism of inflation in this model is 
as follows. In the very early universe before inflation, electromagnetic quantum fluctuations 
are generated due to the breaking of the conformal invariance of the electromagnetic field 
and they act as a source for inflation. Furthermore, also during inflation electromagnetic 
quantum fluctuations are newly generated and the scale is stretched due to inflation, so 
that the scale can be larger than the Hubble horizon at that time, and they lead to the 
large-scale magnetic fields observed in galaxies and clusters of galaxies. This idea is based 
on the assumption that a given mode is excited quantum mechanically while it is subhorizon 
sized and then as it crosses outside the horizon "freezes in" as a classical fluctuation [23]. 
Furthermore, we demonstrate that both inflation and the late-time acceleration of the uni- 
verse can be realized in a modified Maxwell- _F(i?) gravity proposed in Ref. [13] which is 
consistent with solar-system tests and cosmological bounds and free of instabilities. We also 
consider classically equivalent form of non-minimal Maxwell- gravity, and propose the 
origin of the non-minimal gravitational coupling function based on renormalization-group 
considerations. 

This paper is organized as follows. In Sec. II we consider a non-minimal gravitational 
coupling of the electromagnetic field in general relativity. First, we describe our model and 
derive equations of motion from it. Next, we consider the evolution of the large-scale electric 
and magnetic fields. Furthermore, we analyze the gravitational field equation, and then show 
that power-law inflation can be realized. In Sec. Ill we consider a non-minimal gravitational 
coupling of the electromagnetic field in a modified gravitational theory proposed in Ref. [13]. 
We show that in this theory both inflation and the late-time acceleration of the universe 
can be realized. In Sec. IV we consider classically equivalent form of non-minimal Maxwell- 
F(R) gravity. Finally, some summaries are given in Sec. V. In Appendix, we propose the 
origin of the non-minimal gravitational coupling function based on renormalization-group 
considerations. 

We use units in which fee = c = H — 1 and denote the gravitational constant 8irG by k 2 , 
so that k 2 = 8tt/M P i 2 , where M P i = G~ 1/2 = 1.2 x 10 19 GeV is the Planck mass. Moreover, 
in terms of electromagnet ism we adopt Heaviside-Lorentz units. 

II. INFLATION IN GENERAL RELATIVITY 

In this section, we consider a non-minimal gravitational coupling of the electromagnetic 
field in general relativity. 

A. Model 

We consider the following model action: 



In Ref. [28], gravitational-electromagnetic inflation from a 5-dimensional vacuum state has been consid- 






(2.2) 



crcd. 
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C EM = -\l(R)F llu F^\ (2.3) 

I(R) = 1 + f(R) , (2.4) 

where g is the determinant of the metric tensor g^, R is the scalar curvature arising from 
the spacetime metric tensor g^, and £eh is the Einstein-Hilbert action. Moreover, F^ v = 
d^Ay — d u An is the electromagnetic field-strength tensor. Here, is the U(l) gauge field. 
Furthermore, f(R) is an arbitrary function of R. 

The field equations can be derived by taking variations of the action Eq. (2.1) with respect 
to the metric g^ u and the U(l) gauge field A^ as follows: 

itW - \%uR = , (2.5) 



with 



and 



+ ^f{R)F aP F^R llv + g, u U [f'(R)F a(3 F a ^] - V M V, [f(R)F aP F^] J , (2.6) 



~ -.8, {V^I(R)F^) = , (2.7) 



where the prime denotes differentiation with respect to R, V M is the covariant derivative 
operator associated with g^, and □ = g^V^Vj, is the covariant d'Alembertian for a scalar 

field. In addition, R^ v is the Ricci curvature tensor, while T^ M ^ is the contribution to the 
energy-momentum tensor from the electromagnetic field. 

We assume the spatially flat Friedmann-Robertson- Walker (FRW) spacetime with the 
metric 

ds 2 = -dt 2 + a 2 (t)dx 2 = a 2 (r])(-dr] 2 + dx 2 ) , (2.8) 

where a is the scale factor, and 77 is the conformal time. In this spacetime, g^ v = 
diag (— 1, a 2 (t), a 2 (t), a 2 (t)), and the components of R^ and R are given by 

i?oo = -3 (H + # 2 ) , R 0l = , Rij = (h + 3# 2 ) g i3 , R = q(h + 2# 2 ) , (2.9) 

where H = a/a is the Hubble parameter. Here, a dot denotes a time derivative, ' = d/dt. 



B. Evolution of large-scale electric and magnetic fields 

First, we consider the evolution of the U(l) gauge field in this background. Its equation 
of motion in the Coulomb gauge &> Aj(t, x) — and the case of Ao(t, x) = 0, reads 



I\ ■ 1 (3) 

Ai(t, x) + ( H + j \ Ai(t, x) - —A Ai(t, x)=0, (2.10) 
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(3) 

where A = d l di is the flat 3- dimensional Laplacian. 

We shall quantize the U(l) gauge field A^(t,x). It follows from the Lagrangian of the 
electromagnetic field (2.3) that the canonical momenta conjugate to A^(t,x) are given by 



7T = 0, 7Tj = Ia(t)Ai(t, x) . 
We impose the canonical commutation relation between Ai(t,x) and 7Vj(t,x) 



[Ai(t,x),-Kj(t,y) \ = i J 



d 3 k ik . (x _ y) ( , _ hki 

(27T) 3 V 3 & 



(2.11) 



:2.12) 



where k is comoving wave number and k = \k\. From this relation, we obtain the expression 
for Ai(t, x) as 

—3^ Yl U(fc^)e i (fc^)^,t)e lfe - + 6 t (fc,a)e*(fc,a)A*(A ; ,t)e- ife - ,(2.13) 

( 27r ) a=l,2 

where 6j(fc, a) (a = 1, 2) are the two orthonormal transverse polarization vectors, and b(k, a) 
and b^(k,a) are the annihilation and creation operators which satisfy 



b(k, a), [k,a 



b(k,a),b(k,a) = & (k,a),bUk,d) =0.(2.14) 



It follows from Eq. (2.10) that the mode function A(k,t) satisfies the equation 

A(k,t) + \H + j) A(k,t) + ^A(k,t) = 0, 



(2.15) 



and that the normalization condition for A(k, t) reads 



A(k,t)A*(k,t) - A(k,t)A*(k,t) 



la 



Replacing the independent variable t by rj, we find that Eq. (2.15) becomes 



d 2 A(k,r]) 1 dI(r))dA(k,r)) 



Or] 2 



I (if) drj drj 



+ k 2 A(k,r]) =0. 



(2.16) 



(2.17) 



We are not able to obtain the exact solution of Eq. (2.17) for the case in which I is given 
by a general function of rj. In fact, however, we can obtain an approximate solution with 
sufficient accuracy by using the Wentzel-Kramers-Brillouin (WKB) approximation on sub- 
horizon scales and the long-wavelength approximation on superhorizon scales, and matching 
these solutions at the horizon crossing [29, 30]. 

In the exact de Sitter background, we find —krj = k/(aH). Moreover, at the horizon 
crossing, H = k/a is satisfied, and hence — kr]k = 1 is satisfied. Here, rjk is the conformal 
time at the horizon-crossing. The subhorizon (superhorizon) scale corresponds to the region 
A; 1 77 1 ^> 1 (k\rj\ 1). This is expected to be also a sufficiently good definition for the 
horizon crossing for power-law inflation a oc t p , where p> 1, which is almost equivalent to 
exponential inflation because in this case, — ki) = \pj (p — 1)] kj (aH) w kj (aH). 
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The WKB subhorizon solution is given by 

1 



r 1/2 (v)e 



■ikr) 



(2.18) 



where we have assumed that the vacuum in the short-wavelength limit is the standard 
Minkowski vacuum. 

On the other hand, the solution on superhorizon scales, A ont (k, rj), can be obtained by 
using the long- wavelength expansion in terms of k 2 and matching this solution with the WKB 
subhorizon solution in Eq. (2.18) at the horizon crossing. The lowest order approximate 
solution of 74 out (A;,?7) is given by [29] 



4>ut(M) = C(k) + D(k) 



I(rj) 



(2.19) 



where 



C(k) 
D(k) 



i- 1/2 (v) 



1 dl{rj) 

2 drj 



+ ikl(r)) 



) r 


e -ikr) 


I J, 1(1,) 





V=Vk 



l —rV\rf) ( l^P- + ikl(r))) 



2k 



2 drj 



(2.20) 
(2.21) 



r)=Vk 



Neglecting the decaying mode solution, from Eqs. (2.19) and (2.20) we find that \A(k,r))\ 2 
at late times is given by 



\A(k, v )\ 2 =\C(k)f 



2kI( Vk 



1 1 dljrjk) 
2kI( Vk ) drj 



+ i 



Vk 



[2.22) 



where rjf is the conformal time at the end of inflation. 
The proper electric and magnetic fields are given by 

^ proper (t,cc) = a^Eifax) = -a^Aifax) , 
Br opw (t,x) = ar x Bi(t,x) = a~ 2 e l3k d 3 A k (t,x) 



^2.23) 
^2.24) 



where Ei(t,x) and B{(t,x) are the comoving electric and magnetic fields, and is the 
totally antisymmetric tensor (ei23 = 1). 
Using Eqs. (2.19) and (2.23), we find 



\E^(k, V )\ 2 = 2- 



dA(k,r)) 



v \m\ 2 



(2.25) 



where the factor 2 comes from the two polarization degrees of freedom. Multiplying 
\E proper (k, r))\ 2 in Eq. (2.25) by the phase-space density, ink 3 / (27r) 3 , we obtain the amplitude 
of the proper electric fields in the position space 



Airk 3 

|£Propc r(jM) |2 = __|£ Pr opc r(M)| 2 



\D(k)\ 2 k 4 



1 



ir 2 k a 4 |/(r/)|' 



(2.26) 



6 



on a comoving scale L = 2it/k. Furthermore, the energy density of the large-scale electric 
fields in the position space is given by 



M£) , )= l| B ™- (i , l , ) |V(,) = M!^. (2.27) 
Similarly, using Eqs. (2.19) and (2.24), we find 

\B^(k,rj)\ 2 = 2^\A(k lV )\ 2 = 2^\C(k)\\ (2.28) 

where the factor 2 comes from the two polarization degrees of freedom. Multiplying 
|£> proper (fc, 77) | 2 in Eq. (2.28) by the phase-space density, Airk 3 / {2tt) 3 , we obtain the amplitude 
of the proper magnetic fields in the position space 

| B p r o Per(jM)|2 = ^\B^(k,rj)\ 2 = (2-29) 

on a comoving scale L = 2it/k. Furthermore, the energy density of the large-scale magnetic 
fields in the position space is given by 



Pb(L, V ) = ^\B pTopeT (L,r])\ 2 I(r)) = W^^Ifr) . (2.30) 



Here we note the following point. As an example, if / is given by the following form: 
I(rj) = I s (f]/rj s )~ a , where r] s is some fiducial time during inflation, J s is the value of I(rj) at 
77 = ?7 S , and a is a constant, from Eq. (2.20) we find k\C{k)\ 2 = C / [2J(%)], where C is a 
constant of order unity [29, 30]. Hence, using this relation and Eq. (2.30), we find 

c fk\ 4 i(n) 



PB (L,Tl) = —p[-\ j^r. (2.31) 

{2ny \aj I{Vk) 

C. Power-law inflation 

The (fJi,v) = (0,0) component and the trace part of the (//, v) = component of 

Eq. (2.5), where % and j run from 1 to 3, read 



H 2 + J l 



y {/(/?) (g a "F op F 0a - l -g m F aP F a ^j 

f(R) (H + H 2 ^j + 6f"(R)H [H + AHH\ 



3 r 

+ 2 



3 1 1 (3) 1 

+ -f(R)H {FapF*)* - g/W^A (F aP F^) j , (2.32) 
J x = -F(R) - F'(R) (H + H 2 ) , (2.33) 
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and 



2H + 3H 2 + Jo 



K f 1 
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: I{R)F aP F af) 



f(R) [H + 3H' 



+ 6f"(R) \ H + 7HH + 4/j 2 + 12H 2 H) + 3Qf"(R) \ H + 4M 
+ 3 [/'(fl)tf + 4f"(R) (H + 4m)] (F aP F a Py + /'(i?) (F Q/3 F Q/3 )' 



F a pF a P 



[2.34) 



J 2 = -F(R) - F'(R) (H + 3H 2 )+ QF"(R) 



H + 4:1 H 



1 HHJ 



where 



+ 36F'"(R) [H + 4HH 



g a/3 F op F 0a - \gwF aP F^ = i (\F 



proper 



F a pF al3 = 2(\Br° pcr (t,x) 



\Br opeT (t,x)\ 2 ) 
\Er oper {t,x)\ 2 ) 



(2.35) 

(2.36) 
(2.37) 



respectively. Here, a large dot in terms of F a pF a ^ denotes a time derivative, = 
d (F aP F al3 ) /dt. Moreover, Ji and J 2 are correction terms in a modified gravitational theory 
described by the action in Eq. (3.1) in the next section. Hence, because in this section we 
consider general relativity, i.e., the case F(R) = in the action in Eq. (3.2), here both J\ and 
J2 are zero. In deriving Eqs. (2.32) and (2.34), we have used equations in (2.9). Moreover, in 
deriving Eqs. (2.36) and (2.37), we have used Eqs. (2.23) and (2.24). Furthermore, applying 
Eqs. (2.26) and (2.29) to |/^ proper (t, x)\ 2 and \Bi pTopeI (t, x)\ 2 , respectively, we find 



[F a0 F^ 



-H\B pvopci (L,ri)\< 



+ 



H + 3 



f\R) 
1 + f(R) 



H 



\E propex (L,ri)\' 



(2.38) 



Here we consider the case in which magnetic fields are mainly generated rather than 
electric fields because we are interested in the generation of large-scale magnetic fields. It 
follows from Eqs. (2.27) and (2.30) that this situation is realized if / increases rapidly in time 
during inflation [30]. (Hence, from this point we neglect terms in electric fields.) Moreover, 

(3) 

we consider the case in which A {F al3 F a f) is very small because it corresponds to the second 
order spatial derivative of the quadratic quantity of electromagnetic quantum fluctuations, 
so that it can be neglected. In this case, using Eqs. (2.29) and (2.38), we find that Eqs. (2.32) 
and (2.34) are reduced to 



H 



K 



^I(R) - f'(R) [H + hH 2 )+ 6f"(R)H ( H + 4.HH 



k\C(k)\ 2 k 4 



Ti 



[2.39) 
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and 



2H + 3H 2 = k 2 



I(R) + f\R) (-5H + H 2 ) + 6f"(R) [H — HH + AH 2 - 20H 2 H 



36/'" (R) [H + AHH 



k\C(k)\ 2 k A 
V 2 a 4 



respectively. Eliminating I(R) from Eqs. (2.39) and (2.40), we obtain 



(2.40) 



H + H 



f'(R) (-2H + 3H 2 ) +3f"(R) ( H - 2HH + AH 2 - 2AH 2 H 



18 f"'(R) (H + AHH 



k\C(k)\ 2 k 4 
7 2 ^ 



(2.41) 



Here we consider the case in which f(R) is given by the following form: 

ci (R/m 2 ) n 



f(R) = fns(R) = 



c 2 {R/m 2 ) n + 1 



which satisfies the conditions: 



lim fns{R) — — — const 

R— >oo C2 

gm/Hs(fl) = 0. 

/i — 



(2.42) 

(2.43) 
(2.44) 



Here, ci and C2 are dimensionless constants, n is a positive constant, and m denotes a mass 
scale. This form, fns(R), has been proposed by Hu and Sawicki [8]. The second condition 
(2.44) means that there could exist a flat spacetime solution. Hence, because in the late 
time universe the value of the scalar curvature becomes zero, the electromagnetic coupling 
/ becomes unity, so that the ordinary Maxwell theory can be naturally recovered. 

In order to show that power-law inflation can be realized, we consider the case in which 
the scale factor is given by a(t) = a {t/t) p , where t is some fiducial time during inflation, a 
is the value of a{t) at t = t , and p is a positive constant. In this case, H = p/t, H = —p/t 2 , 
H = 2p/t 3 , and H = —Qp/t 4 . Moreover, it follows from the fourth equation in (2.9) that 
R = 6p(2p — l)/t 2 . At the inflationary stage, because R/m 2 3> 1, we are able to use the 
following approximate relations: 



fns(R) 

fus(R) 
fus(R) 



Cl 

c 2 



1 

c 2 



R 



nci 1 / R 



cz 2 m 2 \m 2 



-(n+l) 



n(n + l)ci 1 



c 2 < 



m 4 V m 2 



R 



-(n+2) 



n(n + l)(n + 2)ci 1 / R 



rrr \ m 



-(n+3) 



(2.45) 

(2.46) 
(2.47) 
(2.48) 
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Substituting the above relations in terms of a, H, and R and Eqs. (2.46)-(2.48) into 
Eq. (2.41), we find 



n + 1 



(2.49) 



P = 



2 



a 




1/4 



f^k\C{k)\ 2 k 4 K 2 m 



2n 



(2.50) 



Hence, if n 3> 1, p becomes much larger than unity, so that power-law inflation can be 
realized. Consequently, it follows form this result that the electromagnetic field with a 
non-minimal gravitational coupling in Eq. (2.3) can be a source of inflation. 

Here we state the two following points. In this paper we consider only the case in 
which the values of the terms proportional to f'(R), f"{R) and f"'(R) in the right-hand 
side of Eqs. (2.39) and (2.40) are dominant to the value of the term proportional to I(R). 
Among the terms proportional to f'{R), f"{R) and f"'(R), the term proportional to f'(R) is 
dominant, and its value is order f'(R)H 2 m n {c\/c2 2 ) (H 2 /m 2 ) (R/m 2 ) < - n+1 ' ) ) which follows 
from Eq. (2.46). Here, it follows from H = p/t and R = Qp(2p— l)/t 2 that R is order 10H 2 . 
The condition that the term proportional to f'(R) is dominant in the source term would be 
I(R)/ [f'(R)H 2 ] ~ 10c2(R/ m 2 ) n / n <C 1. This would require extremely small C2 because at 
the inflationary stage R/m 2 3> 1 and n> 1. In such a case, the value of the right-hand side 
of Eq. (2.41), which is order n 2 f'(R)H 2 pB/I (this estimation is derived by using Eq. (2.30)), 
can be order H 2 . Consequently, the right-hand side of Eq. (2.41) can balance with the left- 
hand side of Eq. (2.41), and hence Eq. (2.41) can be satisfied without contradiction to the 
result, i.e., power-law inflation in which p is much larger than unity can be realized. The 
reason why we consider the case in which the term proportional to I(R) on the right-hand 
side of Eqs. (2.39) and (2.40) is so small in comparison with the term proportional to f'(R) 
that it can be neglected is as follows: If the opposite case, namely, the term proportional to 
I{R) is dominant to the term proportional to f'(R), Eqs. (2.39) and (2.40) are approximately 
written as H 2 « k 2 pb/6 and 2H + 3H 2 ps k 2 pb/6, respectively. Thus, in this case it follows 
from Eqs. (2.39) and (2.40) that H 2 and 2H + 3H 2 are the same order and their difference, 
2H + 2H 2 , must be much smaller than H 2 . In fact, Eq. (2.41) implies that H + H 2 balances 

with much smaller quantity than k 2 Pb- Now, yH + H 2 ^j / H 2 = (p — l)/p and hence p must 

be very close to unity. Consequently, in this case power-law inflation cannot be realized. 

Furthermore, when we consider the non-minimal electromagnetic theory described by 
Eq. (2.3), in the very early universe before the beginning of inflation electromagnetic quan- 
tum fluctuations can be generated due to the breaking of the conformal invariance of the 
electromagnetic field through its non-minimal gravitational coupling, This is because it is 
considered that in the very early universe before inflation (e.g., the grand unified theory 
(GUT) scale), there can exist quantum fluctuations of all physical quantities, as the quan- 
tum fluctuations of the inflaton field in the chaotic inflation scenario [31]. On the other 
hand, the non-minimal coupling between the electromagnetic field and the scalar curvature 
function f(R) is purely classical. Furthermore, as explained above, in this paper we consider 
the case in which the term proportional to f'(R) in the right-hand side of Eqs. (2.39) and 
(2.40) is dominant to the term proportional to I(R)- Hence, power-law inflation can be 
realized due to not the term proportional to I{R), namely, the energy density of large-scale 
magnetic fields, but the term proportional to f'{R), namely, a non-minimal electromagnetic 
coupling. Consequently, in this model we consider that inflation can be realized due to not 
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purely quantum effects but semi-classical effects. 

Finally, we note the following three points. In the present model, large-scale magnetic 
fields can be generated due to the breaking of the conformal invariance of the electromagnetic 
field through a coupling with the scalar curvature, I{R)F fJLV F^ v , as is shown in the preceding 
subsection. If there does not exist such a coupling, i.e., f(R) = and hence J = 1, in which 
the ordinary Maxwell theory is realized, electromagnetic quantum fluctuations cannot be 
generated in the FRW spacetime because this background spacetime is conformally flat. 
This result is also realized in the case of dilaton electromagnetism [32-36], in which the 
Lagrangian of the electromagnetic field is given by I{^)F^ V F^ U with /($) = e XK ^ [35], 
where $ is the dilaton field and A is a dimensionless constant. (It is also realized in other 
scalar-field electromagnetism [37-39].) 

Moreover, Bertolami and Paramos have recently considered constraints on a non-minimal 
gravitational coupling of matter, namely, for the present model, f(R) in Eq. (2.4), from the 
observational data of the central temperature of the Sun [20] . They have studied the effect 
of a non-minimal gravitational coupling of matter on the hydrostatic equilibrium of the 
spherically symmetric system with a polytropic equation of state approximately describing 
the Sun with sufficient accuracy, assuming a perturbative regime to the usual Tolman- 
Oppenheimer-Volkoff (TOV) equation of hydrostatic equilibrium and taking into account the 
validity of the Newtonian regime in a theory with a non-minimal gravitational coupling of 
matter. According to them, there exists no strong constraints on a non-minimal gravitational 
coupling of matter obtained from the comparison of the predictions of the theoretical models 
and the current observational sensitivity to the central temperature of the Sun except for the 
relation of the perturbative approach, <C 1. It follows from the relation \ f{R)\ <C 1 

that for the case f(R) = fa s (R) in Eq. (2.42) with m = m e = 0.51 lMeV [22], where m e 
is the electron mass, using the maximum value of the central mass density of the Sun, 
p c = 1.62 x 10 2 gcm -3 , and the expression of the scalar curvature in the Newtonian regime, 
R m —87rGp c , we find that the constraint on fns(R) is given by 



I/hs(J2)| 



ci (SirGpc/miy 



c 2 (-8irG Pc /m 2 e ) n + 1 



Ci (-4.51 x 10 



c 2 (-4.51 x 10" 46 ; 



<1. (2.51) 



Furthermore, the existence of the non-minimal gravitational coupling of the electromag- 
netic field f(R) in Eq. (2.4) changes the value of the fine structure constant, i.e., the strength 
of the electromagnetic coupling. Hence, the deviation of the non-minimal electromagnetism 
from the usual Maxwell theory can be constrained from the observations of radio and optical 
quasar absorption lines [40], those of the anisotropy of the cosmic microwave background 
(CMB) radiation [41, 42], those of the absorption of CMB radiation at 21 cm hyperfine 
transition of the neutral atomic hydrogen [43], and big bang nucleosynthesis (BBN) [44, 45] 
as well as solar-system experiments [46] (for a recent review, see [47]). 



III. INFLATION AND LATE-TIME COSMIC ACCELERATION IN MODIFIED 
GRAVITY 

In this section, we consider a non-minimal gravitational coupling of the electromagnetic 
field in a modified gravitational theory proposed in Ref. [13]. 
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A. Inflation 



We consider the following model action: 

Smg = J d A x^g [ Cmg + £em ] , (3.1) 

C MG ^^[R+F(R)] , (3-2) 

where F(R) is an arbitrary function of R. Here, £em is given by Eq. (2.3). We note 
that F(R) is the modified part of gravity, and hence F(R) is completely different from the 
non- minimal gravitational coupling of the electromagnetic field f(R) in Eq. (2.4). 

Taking variations of the action Eq. (3.1) with respect to the metric g^, we find that the 
field equation of modified gravity is given by [13] 

[1 + F'(R)} R^ - l -g^ v [R + F(R)} + g^UF\R) - W^J V F'{R) = k 2 T^ . (3.3) 

The (/i, v) = (0,0) component and the trace part of the (//, v) = component of 

Eq. (3.3), where % and j run from 1 to 3, are given by Eqs. (2.32) and (2.34), respectively. 
Similarly to the preceding section, we here consider the case in which terms in electric fields 

(3) 

and A {Fa0F al3 } are negligible. In this case, eliminating I(R) from Eqs. (2.32) and (2.34), 
we obtain 

H + H 2 + !^F(R) - F'(R)H 2 + 3F"(R) H + 4 (# 2 + HH) + 18F"'(R) (h + AHH^^ 

f'(R) (-2H + 3H 2 ^j + 3f"(R) (h - 2HH + AH 2 - 2AH 2 H s j 

k\C(k)\ 2 k A 



K 2 



+ 18 f"'{R) (h + AHH^ 



it 2 a 4 



(3.4) 



Here we consider the case in which F(R) is given by 



F ( R) = —m 2 !fbMt!iw«r (3JS) 

cz + U {[(R/liP) - (R /AP)f +1 + (R /APf +1 ] 

which satisfies the following conditions: lim^oo F(R) = —M 2 jc^ = const, lim^o F(R) = 
0. Here, C3 and C4 are dimensionless constants, I is a positive integer, and M denotes a mass 
scale. We consider that in the limit R — > 00, i.e., at the very early stage of the universe, F(R) 
becomes an effective cosmological constant, and that at the present time F(R) becomes a 
small constant, namely, 

lim F(R) = -M 2 — = -2Ai , (3.6) 

, s o (RJM 2 ) 21+1 , . 

F(Ro) = —M =rr7- = -2R* , 3.7 

c 3 + c 4 (Ro M* f +1 
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where Ai (p> H 2 ^ is an effective cosmological constant in the very early universe and 
Rq (~ i^o 2 ) is current curvature. Here, H is the Hubble constant at the present time [48]: 
H = 100/i km s" 1 Mpc" 1 = 2.1h x 10" 42 GeV « 1.5 x l(T 33 eV, where we have used 
h = 0.70 [49]. From Eqs. (3.6) and (3.7), we find 



c 3 = 



Rq 
M 2 



21 



^ Rq 



Rq 

M 2 



2; 



c 4 



1M 2 



(3.8) 
(3.9) 



where the last approximate equality in Eq. (3.8) follows from (i?o/A;) <C 1. 

Furthermore, we consider the case in which f(R) is given by the following form: 



f(R) = /no(-R) 



f 2\2<?+l 



[(i?/M 2 ) - (R /M 2 )} 2q+1 + (i? /M : , 
c 5 + c 6 {[(R/M*) - (iyM 2 )] 2<z+1 + (Ro/M^ 1 } 



(3.10) 



which satisfies the following conditions: liniR^oo /no(-R) — l/c6 — const, lim^ /no(-R) = 0. 
Here, C5 and C6 are dimensionless constants, and q is a positive integer. The form of F(R) 
in Eq. (3.5) and f^o(R) in Eq. (3.10) is taken from Ref. [13]. This form corresponds to the 
extension of the form of fns(R) i n Eq. (2.42). It has been shown in Ref. [13] that modified 
gravitational theories described by the action (3.2) with F(R) in Eq. (3.5) successfully pass 
the solar-system tests as well as cosmological bounds and they are free of instabilities. 

At the inflationary stage, because R/M 2 3> 1 and R/M 2 3> Ro/M 2 , we are able to use 
the following approximate relations: 



F(R) w —M 2 — 



("I 



1-^ 

c 4 



R_ 

M 2 



-{21+1) 



(3.11) 



and 



ho(R) « - 
c 6 



_ / p \ -(29+1)' 

c 6 V^ 2 



(3.12) 



At the very early stage of the universe, because R — > cxd, it follows from Eq. (3.6) and 
the condition, lim^oo /no(-R) = V c 6 = const, that Eq. (3.4) are reduced to 



r2 Ai 



H + H*=—. 

3 



(3.13) 



From this equation, we obtain 



a(t) oc exp ( \j -^t ) . 



(3.14) 



Hence exponential inflation can be realized. Thus, we see that the terms in F(R) on the 
left-hand side of Eq. (3.4), i.e., the part of the braces {}, can be a source of inflation, in 
addition to the terms in f(R) on the right-hand side of Eq. (3.4). In fact, if there do not 
exist any terms in F(R), in which the theory is general relativity, or the contribution of the 
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terms in F(R) to inflation is much smaller than those in f(R), Eq. (3.4) is equivalent to 
Eq. (2.41). In such a case, similarly to the consideration in Sec. II C, substituting a(t) oc t^, 
where ft is a positive constant, the approximate expressions of f^ (R), f^ Q (R) and f^' Q (R) 
derived from Eq. (3.12) into Eq. (3.4), we find ft = q + 1. Hence, if q ^> 1, p becomes 
much larger than unity, so that power-law inflation can be realized. Consequently, in the 
present model there exist two sources of inflation, one from the modified part of gravity, 
F(R), and the other from the non- minimal gravitational coupling of the electromagnetic 
field, f(R)- We here note that even if the value of A; is so small that the modification 
of gravity cannot contribute to inflation, inflation can be realized due to the non-minimal 
gravitational coupling of the electromagnetic field, namely, the change of the value of f(R) 
in terms of R, and the generation of magnetic fields. This is an important feature of the 
present model. 



B. Late-time cosmic acceleration 

Next, we consider the late-time acceleration of the universe. As shown above, at the early 
stage of the universe, at which the curvature is very large, inflation can be realized due to the 
terms in F(R) and/or those in f{R)- As curvature becomes small, the contribution of these 
terms to inflation becomes small, namely, the values of these terms in Eq. (3.4) become small, 
and then inflation ends. After inflation, radiation becomes dominant, and subsequently 
matter becomes dominant. When the energy density of matter becomes small and the value 
of curvature becomes Ro, there appears the small effective cosmological constant at the 
present time as seen in Eq. (3.7). Hence, the current cosmic acceleration can be realized. It 
has been shown in Ref. [13] that both inflation and the late-time acceleration of the universe 
can be realized in modified gravitational theories described by the action (3.2) with F(R) in 
Eq. (3.5) for the case without the non-minimal gravitational coupling of the electromagnetic 
field f(R) in Eq. (2.4). In this subsection, we confirm that also in this theory with the non- 
minimal electromagnetic coupling f(R), the late-time acceleration of the universe can be 
realized. (Incidentally, it has been shown in Ref. [21] that in this theory with a non-minimal 
coupling with the kinetic term of a massless scalar field, the late-time acceleration of the 
universe can be realized.) 

In the limit R — *> Rq, i.e., the present time, because R/M 2 — R /M 2 C 1, we are able to 
use the following approximate relations: 



F(R) 



-M 



c 3 



c 3 + c 4 (R /M 2 



21+1 




Rq 

M 2 



21+1 



c 3 + c 4 (Rp/M 
c 3 



2,21+1 



Rq 

M 2 



21+1 



(3.15) 



and 



fao(R) 



c 5 



[c + cURo/M 2 ) 2 ^ 1 } 2 



x 



(R_ Ro_\ 
\M 2 M 2 J 



2q+l 



c 5 + c 6 (Ro/M 2 ) 2q+1 



c 5 



M 
M 2 J 



2g+l 



(3.16) 



14 



From Eqs. (3.4), (3.15), the approximate expressions of F'(R), F"(R) and F'"(R) derived 
from Eq. (3.15), and the approximate expressions of f^o(R), f^o(R) and f^' Q (R) derived 
from Eq. (3.16), we see that if q > I, /no(-R) becomes constant more rapidly than F(R) in 
the limit R — > Rq. As a result, the electrodynamics looks as purely minimal theory at the 
current universe. For such a case, in the limit R — > Rq, Eqs. (3.4) are reduced to 



From this equation, we obtain 



H + H 2 = ^-. (3.17) 



so that 

I = f>°- ^ 

Thus, the late-time acceleration of the universe can be realized. 

Finally, we note the following point: In this model, even if the value of Rq is so small that 
the modification of gravity cannot contribute to the late-time acceleration of the universe, the 
late-time acceleration can be realized due to due to the non-minimal gravitational coupling 
of the electromagnetic field and the generation of magnetic fields. This is also an important 
feature of the present model. 



IV. CLASSICALLY EQUIVALENT FORM OF NON-MINIMAL MAXWELL-F(^) 
GRAVITY 

In this section, we consider classically equivalent form of non-minimal Maxwell- _F(i?) 
gravity. 

The action (3.1) can be rewritten by using auxiliary fields. Introducing two scalar fields 
£ and £, we can rewrite the action (3.1) to the following form [15, 21]: 

S = J dV=0 | JL [c + F(C)} + l(()C M + £ (R - C) } , (4-1) 

C M = ~F^\ (4.2) 

where Cm is the Lagrangian describing the ordinary Maxwell theory. The form in Eq. (4.1) 
is reduced to the original form in Eq. (3.1) by using the equation ( = R, which is derived by 
taking variation of the action (4.1) with respect to one auxiliary field £. Moreover, taking 
variation of the form in Eq. (4.1) with respect to the other auxiliary field £, we find 

e=^[i+nc)]+/'(c)£ M , (4.3) 
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where the prime denotes differentiation with respect to £. Substituting Eq. (4.3) into 
Eq. (4.1) and eliminating £ from Eq. (4.1), we find 



S = J dV=${ JL [i + f\C)] R + [/(C) + /'(C) (R - 0] £m 



+ ^ [*X0 - nOC] 



We make the following conformal transformation of the action given by Eq. (4.4): 



with 



e* = l + F'(C) 



(4.4) 



(4.5) 



(4.6) 



where tp is a scalar field. Here, the hat denotes quantities in a new conformal frame in 
which the term in the coupling between F'(Q and the scalar curvature in the first term on 
the right-hand side of Eq. (4.4) disappears. Consequently, the action in the new conformal 
frame is given by [50] 



d — g 



± (r - \r v d^d vV 



+ ( e" 2 ^ {I [CM] - I' [CM] CM} + e*I' [CM] [R + aOp- -grd&dv* 



c 



M 



^ e - 2 n^[CM]-K-i)CM} 



(4.7) 



where 



(4.8) 



and g is the determinant of g^ u . In deriving Eq. (4.7), we have used Eq. (4.6). Moreover, 
CM in Eq. (4.7) is obtained by solving Eq. (4.6) with respect to C as C = CM- Hence, the 
action in the new conformal frame (4.7) includes the Brans-Dicke type scalar field <p [51]. 
From the term proportional to Cm on the right-hand side of Eq. (4.7), we see that the form 
of the Lagrangian in terms of the electromagnetic field in Eq. (4.7) is close to that of the 
electromagnetic field with the coupling to the dilaton, which has been explained in Sec. II C. 
In other words, the Lagrangian of non-minimal Maxwell- _F(i?) gravity is qualitatively similar 
to Lagrangian describing dilaton electromagnetism (except for the term in the coupling 
between the scalar curvature and the electromagnetic field). As explained in Refs. [52, 53], 
however, this fact does not mean the physical equivalence between them. 
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V. CONCLUSION 



In the present paper, we have considered inflation and the late-time acceleration in the 
expansion of the universe in non-minimal electromagnetism, in which the electromagnetic 
field couples to the scalar curvature function. As a result, we have shown that power-law in- 
flation can be realized due to the non-minimal gravitational coupling of the electromagnetic 
field, and that large-scale magnetic fields can be generated due to the breaking of the confor- 
mal invariance of the electromagnetic field through its non-minimal gravitational coupling. 
Furthermore, we have demonstrated that both inflation and the late-time acceleration of the 
universe can be realized in a modified Maxwell-F(i?) gravity proposed in Ref. [13] which 
is consistent with solar-system tests and cosmological bounds and free of instabilities. We 
have also considered classically equivalent form of non- minimal Maxwell-F(i?) gravity. 

Finally, we make a remark about the observational deviation of a non-minimal electro- 
magnetic theory from the ordinary Maxwell theory. It follows from the fourth equation in 
Eq. (2.9) that in exponential inflation the scalar curvature is proportional to the square of 
the Hubble parameter. Moreover, it is known that the root-mean-square (rms) amplitude 
of curvature perturbations is also proportional to the square of the Hubble parameter. In a 
non-minimal electromagnetic theory, because magnetic fields couple to the scalar curvature, 
there can exist the cross correlations between magnetic fields and curvature perturbations 
through the Hubble parameter. Hence, if the primordial large-scale magnetic fields are de- 
tected [54, 55] by future experiments such as PLANCK [56], SPIDERS (post-PLANCK) [57] 
and Inflation Probe (CMBPol mission) in the Beyond Einstein program of NASA [58] on the 
anisotropy of CMB radiation, and if there exist (do not exist) the cross correlations between 
the primordial large-scale magnetic fields and curvature perturbations, it is observationally 
suggested that at the inflationary stage there should exist a non-minimal gravitational cou- 
pling of the electromagnetic field (the strength and/or the form of non-minimal coupling of 
the electromagnetic field may be observationally restricted). 
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APPENDIX: ASYMPTOTIC FREEDOM VERSUS NON-MINIMAL COUPLING 

It is very interesting that one can generalize the discussion of this work for interacting 
theories: scalar/spinor electrodynamics and non-Abelian gauge theory. As a simple example, 
let us consider the £77(2) gauge theory with the Lagrangian: C = — (1 / 4) G'j lu G a ' JlU , where 
G" is the SU(2) field strength. The effective renormalization-group improved Lagrangian 
for such a theory in matter sector has been found in Ref. [59] for a de Sitter background as 

1 9 2 
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with 



where g(t) is the running SU{2) gauge coupling constant, g is the value of g(t) in the case 
t — 0, and t is a renormalization-group parameter. Note that the running gauge coupling 
constant typically shows asymptotically free behavior: it goes to zero at very high energy. 
For the covariantly constant gauge background with G a ^ u G aflv /2 = H 2 , where H corresponds 
to the magnetic field in the 577(2) gauge theory, it has been proposed in Ref. [59] that t is 
given by 



where fi is a mass parameter. 

It is clear that with the decrease of the energy scale (namely, as the universe expands), 
t is decreasing, as t is very large at the early universe. Taking into account the results of 
this work, one can try to relate the asymptotic freedom in a non-Abelian gauge theory with 
non-minimal Maxwell-modified gravity. In this way, using the proposal of Eq. (2.42) in Sec. 
II for non-minimal f(R) in front of G a v G a ^ v , one gets 

c lW mr nf f (A4) 



c 2 {R/m 2 ) n + 1 12tt 2 

Hence, according to this assumption, at a very large curvature (R/m 2 ^> 1), t ~ 
[127r 2 / (llg 2 )} (ci / c 2 ) , while at the current universe (R — > 0), t — > 0. Thus, asymptotic 
freedom induces the appearance of the non-minimal gravitational gauge coupling in (non-) 
Abelian gauge theories at high energy. 

Generally speaking, such a scenario is universal and it works not only for asymptotically 
free theories. For instance, for scalar QED one can easily write the renormalization-group 
improved effective Lagrangian in curved spacetime. In the matter sector (zero scalar field 
background) it has qualitatively the same form as Eq. (A.l), the only sign of t is different 
in the expression for the running gauge coupling constant. As a result, such an effective 
Lagrangian again induces the non-minimal gravitational coupling of the electromagnetic 
sector. 
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